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Abstract 

Inhomogeneities are introduced in loop quantum cosmology using regular lattice 
states, with a kinematical arena similar to that in homogeneous models considered 
earlier. The framework is intended to encapsulate crucial features of background 
independent quantizations in a setting accessible to explicit calculations of pertur- 
bations on a cosmological background. It is used here only for qualitative insights 
but can be extended with further more detailed input. One can thus see how several 
parameters occuring in homogeneous models appear from an inhomogeneous point of 
view. Their physical roles in several cases then become much clearer, often making 
previously unnatural choices of values look more natural by providing alternative 
physical roles. This also illustrates general properties of symmetry reduction at the 
quantum level and the roles played by inhomogeneities. Moreover, the constructions 
suggest a picture for gravitons and other metric modes as collective excitations in a 
discrete theory, and lead to the possibility of quantum gravity corrections in large 
universes. 



1 Introduction 

Loop quantum gravity PQ 12 El is a candidate for a non-perturbative and background 
independent quantization of general relativity which by now has uncovered several crucial 
features relevant at small length scales. Most importantly, it unambiguously leads to 
discreteness of spatial geometry at the Planck scale [H El [B^ which has played crucial roles 
in physical pictures studied so far, such as cosmological and black hole singularities [HIHlin]- 
These applications, however, also require a significant amount of dynamical aspects which 
remain poorly understood and quite complex in a full setting. Many applications have 
therefore been formulated in symmetric models of loop quantum cosmology piU ITT] , which 
is a reduction preserving crucial properties such as the very discreteness of the full theory 
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but at the same time making calculations explicitly treatable. Results concerning new 
behavior at small scales then indeed follow at least qualitatively as direct consequences of 
such basic properties related to discreteness. 

As homogeneous models in this framework are becoming studied in more and more de- 
tail, initially designed structures have to be refined. In particular, inhomogeneities are the 
crucial ingredient for most remaining open issues. While background independence is con- 
ceptually important and crucial for basic properties exploited in loop applications, it also 
makes more detailed derivations of physical properties complicated. A relation between 
quantum and classical geometry is needed for nearly every aspect considered in quantum 
gravity such as the singularity issue (where a relation to classical geometry is required in 
identifying states where a classical singularity would occur) or corrections to cosmological 
structure formation. While full quantum geometry as it occurs at a fundamental level 
provides many advanced techniques, a direct use usually comes with too much baggage for 
a given application. Additional specializations to pick the right regime are often required. 
Examples are symmetric states in considerations of singularities, or boundary conditions 
for black hole entropy calculations ^21 ^| and in recent attempts to make contact with 
low energy scattering [T3j. One is thus often led to introduce special states corresponding 
to the physical situation at hand, which allows one to make a detailed relation to classi- 
cal properties. These states, in turn, provide useful background structures. As has been 
demonstrated by now in many investigations, characteristic properties of the background 
independent framework remain intact even when such additional structures are introduced 
to extract and focus on physical regimes of interest. Such a narrowing- down is a technical 
step to perform explicit calculations or to find suitable approximation schemes, rather than 
a conceptual part of the underlying theory. Given the complexity as well as incomplete- 
ness of the full framework, a systematic derivation of dynamics of such sectors from the 
full theory is currently difficult, but many qualitative and semi-quantitative aspects can 
be studied in carefully designed models. If such restrictions are not made, immense tech- 
nical difficulties have to be faced. For cosmological perturbation theory for instance the 
averaging problem P3j, which is even difficult in classical gravity, arises while difficulties 
relevant for the singularity problem have been described in [16j. The aim is certainly to 
arrive, at some point, at a complete relation to, if not a derivation from, a fully background 
independent theory. But such a derivation as well as the construction of the full theory 
itself will be much easier if likely properties are already known and have been evaluated in 
simpler situations.^ 

An inclusion of inhomogeneous degrees of freedom in applications is currently under 
construction, and already available for midisuperspace models [TH| CHI 120] and in pertur- 
bative form fSlJ. From these developments it is becoming clear that properties of inho- 
mogeneous degrees of freedom have a bearing also on strictly homogeneous models. While 
the overall findings in the homogeneous case are confirmed, the inhomogeneous analysis 
suggests changes in constructions which one would not necessarily have considered from 

"'^Eddington's observation ". . . that it is indeed helpful in the quest for knowledge if we understand the 
nature of the knowledge we are looking for" JT] also applies in this case. 
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within those models. A relation between homogeneous and inhomogeneous models can in 
particular influence what is seen as natural or unnatural in the constructions, and which 
corrections should be expected to dominate. 

We will discuss here basic properties of quantum variables used in loop models, and 
draw four conclusions regarding parameters of isotropic models when they are viewed as 
arising from inhomogeneous ones. The setting we mainly have in mind is that of pertur- 
bative inhomogeneities and thus refer particularly to the behavior of models when they 
describe larger universes rather than regimes close to classical singularities. Most of our 
observations, however, will be general and can be checked non-perturbatively in midi- 
superspace models, although they are justified by the specific formulation introduced here 
only in regimes of perturbative inhomogeneities. We will also discuss technical details of 
symmetry reductions, suggest a new picture for how metric modes on a background geome- 
try such as gravitons arise effectively through collective quantum excitations, and indicate 
the possibility of quantum corrections to the universe evolution at large scales. 



2 Isotropic loop quantum cosmology 

A classical isotropic model is fully described by the scale factor a{t) as a solution to the 
Friedmann equation once the matter content is specified. In Ashtekar variables fI2\ 1^ on 
which loop quantum gravity is based the scale factor and its time derivative are expressed 
through a densitized triad component p with \p\ = a^, determining an isotropic densitized 
triad E'^ = p6f, and a connection variable c = 70, determining an isotropic connection = 
C(5a [21] (with the Barbero-Immirzi parameter 7 [23, 23] relating the Ashtekar connection 
to extrinsic curvature). The densitized triad encodes spatial geometry as it is related to 
the spatial metric by E^E\ = q"-^ det q. In particular, the spatial volume of a region R 
is V = d^x^/l det E\ which for an isotropic triad reduces io V = Vq\P\^/'^ =: 
where the coordinate volume Vq is absorbed in p. The re-scaled variable p, and similarly 
c = Vq^'^c, is then independent of coordinates in contrast to p or the commonly used scale 
factor a = \/\p\ It does, however, depend on the coordinate size Vq of the region 

integrated over which can be fixed to be the total space in a closed model but needs to 
be chosen in open models where integrations have to be restricted to finite regions. The 
relation of basic variables to metric or connection components then depends on which 
choice is being made here. Nevertheless, Vq will not appear explicitly when one restricts 
oneself to the variables c and p only, for instance the Poisson relations 

, , 87r7G 
^^'^^ = ^ 

become independent of Vq after rescaling: 

, , 87r7G 
{c,p} = • 

Following the loop quantization, these basic objects are represented on a Hilbert space 
with an orthonormal basis of states given as functions of the connection compo- 



3 



nent by (c|/i) = e^^^'^/'^ Exponentials e*'' '^/^, analogous to holonomies of the full theory, 
then act by multiplication, 

7^^\^l) = \^l + ^') (1) 

and the triad component p by a derivative, 

P\^^) = ^i^lM ■ (2) 

This demonstrates how basic properties of the full theory are preserved: Only exponentials 
of c are represented while it is not possible to obtain an operator for an individual compo- 
nent c directly, and the triad operator p has a discrete spectrum since its eigenstates are 
normalizable. From p one directly obtains the volume operator V = 

Using the basic operators one then has to construct more complicated ones relevant 
for dynamics. For matter Hamiltonians, inverse triad components such as | deti?|~^/^ are 
required which cannot simply be obtained by taking an inverse operator of p because this 
inverse does not exist: p has a discrete spectrum containing zero. Nonetheless, one can 
construct operators corresponding to the classical inverse using Poisson identities such as 

m 

1 

R 

for < / < 1 and some fiQ. On the right hand side, only a positive power of p is required 
which is easily available upon quantization. The Poisson bracket will then become a com- 
mutator, resulting in a well-defined and for isotropic models even bounded operator [2H]- 
The main effect is that inverse powers of p occuring in classical equations are replaced by 
regular functions, which can be read off from eigenvalues of the operators of the form 



/io/ 



(3) 



and do not blow up at p = 0. For instance, in a scalar matter Hamiltonian, occurs 
in the kinetic term which is replaced by a regular function d{p) = di{p/p^,) where / is the 
ambiguity parameter above and = |7j7io^p with a second parameter j arising when one 
rewrites expressions in a form mimicking SU(2) holonomies of the full theory, where one can 
choose an irreducible representation for holonomies. Matrix elements in a representation 
j will then be of the form of exponentials used above with exponents between — j/xqC 
and jfioc, the value in corresponding to the fundamental representation j = ^- This 
procedure results in an effective density of the form di{p/p^) = with 

[231301 

Pi{q) = Yl^'-' {{1 + 2)-' {{q + ly^' - \q - 1^') (5) 
- (/ + ly'q ((g + 1)'+^ - sgn(g - l)|g - 1|'+^)) . 
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Similarly, the gravitational part of the Friedmann equation requires reformulations 
because it contains terms linear and quadratic in c while only exponentials of ic, i.e. almost 
periodic functions of c, can be quantized. Thus, the classical expressions will be replaced 
by a function which reduces to the classical one for small c, i.e. small extrinsic curvature 
in a flat model, while giving corrections at larger curvature. This usually also involves free 
parameters as in the simplest case of choosing 

/Iq ^ sin^ fiQC ~ + O(c^). (6) 

In parallel with full constructions p?T| 122] , this can be understood as arising from the trace 
of a holonomy of around a square loop whose edges are integral curves of symmetry 
generators Xj. The square loop holonomy for edges along directions / and J is of the form 

hjhjhj^hj^ = cos^(iyUoc) + 2(1 + 2ejj^TK) sin^(i/ioc) cos^(iyUoc) + (25/j - 1) sin^(i/ioc) 
+4(r/ - rj)(l + 5jj) sin^(i/ioc) cos(|/ioc) 

which, when appearing in a suitable trace, gives rise to a single term 4 sin^(i/ioc) cos^(|/ioc) = 
sin^/xgc |24j. The parameter /iq then determines the ratio of the coordinate length of an 
edge of the loop to V^^^, and it appears in the holonomy through J dtA^ef^Ti which has 

1 /3 

matrix elements /ioK) ^/^ for an isotropic connection. On top of that, one can again 
choose non-fundamental representations for holonomies with the effect of multiplying /iq 
with a spin label j [SHI Olj • 

In the Poisson brackets (jHl) as well as in replacements (jH)) for polynomials in c we use 
holonomies such that it appears most natural to use the same value for /iq here. Its value 
remains undetermined, but the order of magnitude has been estimated by relating it to 
the lowest eigenvalue of the full area operator, corresponding to the area of a square loop 
used in such a basic dynamical move [22] ■ The argument is, however, incomplete because 
the area operator is used to fix a parameter of the reduced constraint although it does not 
appear in the full one as it is formulated now. It seems against the general viewpoint of 
loop quantum cosmology, formulated e.g. in |TT], to invoke the area operator to quantize 
curvature components in a model while it does not occur in the full constraint. Moreover, 
the coordinate area relevant for the loop in holonomies is different from the geometrical 
area quantized by the area operator. 

These construction steps give rise to quantum corrections of different types to classical 
equations, and also to several parameters to choose which are ultimately to be related to the 
full theory. Open questions in this context concern the relative magnitude of correction 
terms, obtained through modifications such as those in ([S]) or dH)) in effective equations 
[SSI ISni ISZl ISHl ISH HO] , with respect to each other, which is important to know for the 
construction of realistic scenarios, and natural ranges of the parameters. While these open 
issues can be avoided to some degree by sufficiently general arguments such as those for 
singularity removal JH] or phenomenology [H], there are additional problems some of 
which have become visible recently: 

• The parameter j labels an SU(2) representation that needs to be chosen for the con- 
struction of inverse triad operators and is sometimes taken to be significantly larger 
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than 1/2 which would be the lowest allowed value corresponding to the fundamental 
representation. This permits one to remain well in a semiclassical regime while still 
having sizable quantum effects. It also justifies, when used only in inverse power cor- 
rections dH), to ignore other correction terms coming from holonomies in ^ because 
they are much smaller in such regimes. It is thus meaningful for technical purposes 
to work with larger values of j and analyze corresponding effects. But it looks "un- 
natural" that the representation label should be large, rather than just the value 
for the fundamental representation, or that different corrections should use different 
j. Moreover, there are arguments related to a type of quantum stabihty ^Ij or the 
physical inner product |^ which may indicate that the fundamental representation 
at least for holonomy corrections ^ is indeed preferred by internal consistency.^ 

• In the construction of the dynamical equations one has to replace ploynomials in 
c by functions depending on exponentials. This is justified when c is small in a 
semiclassical situation, while close to classical singularities there can be strong dif- 
ferences between such functions which are indeed expected when classical gravity 
breaks down. However, the isotropic connection component c can even become large 
in semiclassical regimes where one would not expect strong quantum effects. This 
is the case whenever there is a positive cosmological constant for which the connec- 
tion component behaves as c ~ y/Xa. Thus, even if A is small, a will grow without 
bound and eventually lead to a large c (see Sec. 4.3 of ;,25 and for an early and 
a more recent discussion). Similar effects, depending on initial values, may occur in 
anisotropic models |44j. 

• In a flat model we have to choose a cell of coordinate size Vq, or a compactification 
to a torus of the same size. This is just auxiliary in fiat models while its value is 
fixed in closed ones. However, this value appears in equations through d{p) and in 
physical quantities such as the scale at which a bounce may occur [IHl ESI HI] • This 
is because the rescaled variables p and c, while they are not coordinate dependent, 
depend on the value of Vq entering in the rescaling. 

These effects have been recognized and in some cases ameliorated by amendments 
introduced in isotropic models. For instance, the issue in the presence of a cosmological 
constant can be resolved if one chooses /io not as a constant but as a p-dependent function 
p jlHlllZI. If, e.g., ~ 1/ \/\p\ at large scales, c is always small in semiclassical regimes. 
Moreover, such a choice can remove the Vo-dependence from the bounce scale resulting 
from holonomy corrections jlH]. However, while such a choice is not inconsistent with 
quantization procedures, homogeneous models cannot justify it convincingly. As with fiQ, 
one is invoking the area operator although it does not appear in current versions of the 
full constraint. Moreover, the resulting operator for curvature components becomes triad- 
dependent which is not realized in the full theory either. We will demonstrate in what 

^Those arguments are, however, not as clear-cut as they are sometimes presented. As the authors of 
|34[ I42| mention, their counter-arguments can easily be evaded by working not with a single larger value 
of j but by appropriate sums of operators involving different j. 
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follows that the intuition behind those modifications is nonetheless borne out by lattice 
constructions. 

In fact, isotropic loop quantum cosmology was originally developed for models which 
exhibit typical effects close to classical singularities in order to see how quantum correc- 
tions can lead to better classical behavior. In the meantime, the resulting singularity 
removal mechanism has been extended non-trivially to inhomogeneous models This 
was successful because for this aspect a single spatial patch is already quite typical and 
thus homogeneous models are reliable (while isotropy itself would be too special; see e.g. 
[in]). These models were not intended to be taken too seriously for quantitative aspects 
at larger scales. Although the universe is homogeneous on such scales to a high degree, 
the classical continuum picture requires space to be made of many "atoms" in a discrete 
world. Thus, at the quantum level one must not ignore inhomogeneities especially when the 
universe becomes large. Some of the above difficulties arise precisely from taking isotropic 
models literally at large scales. One can clearly see this for the cosmo logical constant, 
which implies a large integrated trace of extrinsic curvature just because space becomes 
large. This is the case even if the local curvature scale remains small in a semiclassical 
regime. While homogeneous models do not allow much choice in basic variables, which 
will include some form of the total extrinsic curvature, inhomogeneous models are built on 
more local objects. A better situation can then be expected if inhomogeneous models are 
used where basic variables remain small in semiclassical regimes. 

3 Inhomogeneous effects 

The full theory in current form has several different complications not all of which are 
related to inhomogeneities as they arise as modes on a background. 

3.1 Aspects of loop quantum gravity 

Configuration variables are holonomies associated with arbitary curves in space, whose 
values allow one to distinguish between all connections relevant for general relativity. Since 
the connection takes values in su(2), holonomies are elements of the group SU(2). Their 
matrix elements are multiplication operators in the basic representation of loop quantum 
gravity built on states which are gauge invariant products of holonomies associated with 
edges of arbitrary spatial graphs [101 • While these graphs visualize the discrete spatial 
structure of quantum geometry, they can be arbitrarily fine and no explicit short-scale 
cut-off is present in the theory. Graphs can also be of any topology (knotted or linked) and 
can have vertices of high valence. While holonomies implement the connection as basic 
operators, spatial geometry is given by triads, the conjugate momenta. They are quantized 
to flux operators, associated with 2-surfaces in space, which take non-zero contributions 
only from intersection points between the surface and edges of the graph acted on by the 
flux operator. Their spectra are discrete, further implementing the spatial discreteness of 
quantum geometry. 
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These basic operators appear in more complicated ones relevant for dynamics, such as 
matter Hamiltonians for which inverse triad operators are necessary and the gravitational 
Hamiltonian constraint. While such operators can be constructed in well-defined manners 
|321l22j) their actions are highly involved on arbitrary states. Moreover, in particular the 
gravitational constraint usually changes the graph underlying a state it acts on because 
holonomies around closed loops are used to quantize curvature components. Unless such 
loops are lying entirely on the original graph of a state, the action creates new edges and 
new vertices leading to a finer graph. It is a success, however, that such operators can 
be defined in a well-defined manner at all, considering that their analogs in quantum field 
theory on a curved background space-time would suffer from several infinities. This is 
one of the places where background independence and the quantum representation it leads 
to are crucial and imply characteristic properties of resulting theories. As noted before, 
there is no explicit short-scale cut-off in the theory. Finiteness rather results from the fact 
that Hamiltonian operators act on graph states, and each such graph implies a non-local 
representation of the classical fields. 

In addition to the Hamiltonian constraint, there is the diffeomorphism constraint which 
is implemented as finite diffeomorphisms moving graphs in space. Solving it implies that 
invariant results are independent of the embedding of graphs in space. Most or all of this 
freedom needs to be fixed when calculations are to be done using a background geometry 
for a particular physical regime. In the full setting, this can be achieved by picking suitable 
semi-classical states peaked at a geometry corresponding to the desired background (see, 
e.g., inil)- In this way, the full background independent framework is used but special 
states are selected to pick a regime, making use of an additional (background) structure. 
In addition to the classical background, such states depend on typical quantum aspects 
such as the spread of a state or parameters describing the typical fineness of graphs used 
and thus encoding discrete spatial aspects. 

Unfortunately, there are conceptual and technical difficulties in even defining suitable 
semiclassical states, and working with them at this level is highly involved. One of the 
difficulties is, for instance, that holonomies are SU(2)-valued which requires the use of 
lenghty re-coupling identities in doing explicit calculations. On the other hand, results are 
rarely sensitive to all aspects and can often be reproduced with a high level of accuracy 
in simpler constructions (compare, e.g., inS] with ^^). We will now devise a model which 
allows one to formulate a suitably explicit framework for perturbative inhomogeneities 
and other inhomogeneous regimes such as the BKL picture, maintaining the characteristic 
properties of a background independent quantization. 

3.2 Lattice models 

To capture significant inhomogeneous ingredients we introduce regular lattice states with 
spacing io measured in background coordinates assumed to be flat. The motivation is 
to analyze implications of typical aspects of the full theory such as the fact that states 
are based on graphs and properties of Hamiltonian operators. This will allow one to 
read off characteristic modifications to classical expressions which can then be transferred 
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to effective equations. Lattice links are parts of integral curves of a basis of symmetry 
generators of the background, which also provides an orientation of all edges. Although it 
is not necessary for inhomogeneous states, we assume this lattice to be in a cell of finite size 
Vo for comparison with homogeneous models. The number of lattice blocks and vertices is 
then iV = Vo/^Q. Lattices of this type are special cases of fundamental states, but we also 
make use of background structures to be introduced in the construction. This renders the 
usual floating lattice of a background independent fundamental description into a rigid one 
on a background, suitable, e.g., for cosmological considerations. By explicitly introducing 
a background in this manner we bypass more involved reformulations which would define a 
background through relational objects, akin to relational solutions of the problem of time. 

Usually, spin network states are built on graphs, labeled by SU(2) representations 
on edges and contraction matrices in vertices in order to multiply all holonomies, evalu- 
ated in the labeling representation, to a gauge invariant function of the connection. This 
changes when the construction is used on a background, for instance a homogeneous one 
where integral curves of Killing vector fields define the lattice links. For a perturba- 
tive treatment of inhomogeneities classical variables of the form = p\x)6f^6^ and 

— ^A^)^li)^a ~^'^ii^)^a Suitable where the densitized triad and the first part of the 
connection are diagonal. A diagonal densitized triad is obtained by a gauge choice (such as 
scalar modes in longitudinal gauge) making use of the background around which one per- 
turbs. The connection cannot be completely diagonal, however, because this would violate 
the Gauss constraint in an inhomogeneous situation (see also H^j). The non-diagonal 
part containing ipi{x) comes from the spin connection which is determined completely in 
terms of p\x), and possibly spatial derivatives of the shift vector if it is not zero in the 
chosen gauge. The independent degree of freedom ki{x) in the connection thus appears 
in diagonal form (in gauges where the shift vector vanishes this corresponds to extrinsic 
curvature; otherwise there is an additional contribution to the non-diagonal connection 
part). In fact, we have conjugate fields kj{x) and p\x), with 

{~kj{x),p'{y)} = 8njG6-^6{x,y). (7) 

Taking only the diagonal part of A^, holonomies along links are of the form exp{iokjTj) 
in direction I. The diagonal form of basic variables implies that not all the freedom 
of an SU(2) theory has to be dealt with, effectively Abelianizing the framework as in 
homogeneous models jSUllSllISI- Nevertheless, we can consider higher representations of 
SU(2) such that, in the j-representation, a holonomy has matrix elements exp{ij£okj) and 
lower powers. Independent functions on the space of lattice connections are thus labeled 
by integers where v denotes the vertex position and 7 = 1,2,3 the direction of the link 
starting from the vertex using the orientation of symmetry generators as tangent vectors to 
links. Such functions are then of the form exp(i/i/ ,;^o^//2). In homogeneous models, one 
uses real labels /i for representations of the Bohr compactification of the real line, which 
can be seen as arising from a degeneracy between representation labels and edge length 
[2n!- The functions exp(i/xc/2) are then used in isotropic models, which separate the space 
of isotropic connections. Using integer labels on a fixed lattice states does not allow us to 
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separate all classical connections, but this is not required because choosing a lattice means 
that only functions of a certain scale are being probed. 

At this point, we have assumed diagonal metrics also at the inhomogeneous level 
which effectively Abelianizes the theory: instead of SU(2) calculations we simply work 
with U(l), i.e. we can replace complicated SU(2) recoupling relations by multiplications 
of phases.^ This is the reason why we do not need to introduce vertex labels because 
Abelian holonomies are uniquely multiplied to gauge-invariant functions. Although the 
diagonalization is a truncation of the classical theory, it allows full access to perturba- 
tive inhomogeneous degrees of freedom realized classically. This is also the key reason for 
simplifications that makes it possible to do explicit calculations. 

Since this corresponds to a field theory for classical variables p^{x) and ki{x), we have 
many more basic operators which are nevertheless constructed very similarly to those of 
isotropic models. Holonomies along lattice links of the form 

hj^y = exp{i J dtki/2) ^ exp{iiQki{v)/2) = exp{ikj{v)/2) 

ei,v 

absorbing in kf := iokj, and fluxes through lattice sites Sj^^ perpendicular to an edge ej^y 
as in Fig. [T]are given by Fj^y = p\y)d?y = ilp^{v) =: p\v). Compared to isotropic 

1 /3 

models, Vq is thus replaced by £o and the lattice spacing takes over the role of the cell 
size. Now all expressions are dependent on C.q, but this has meaning as the lattice size, 
although it is just the size measured in background coordinates. Nonetheless, £o is not just 
an auxiliary quantity because it is thought of as arising from a fundamental state after 
re-introducing a background rather than being introduced to discretize a theory by hand. 

Poisson relations between holonomies and fluxes introduced here follow, e.g., from a 
mode decomposition 

~ki{x) = ^A;/(K)e*'^-^ 

K 

of the fields kj{x) and p'^{x) on a symmetric background which we still assume to be fiat 
for simplicity. With our box of size Vq, all components of the wave numbers k summed 

— 1/3 

over are kj = 2'kVq n with positive integers n. From 

j hi{x)p\x)&'x = J2~ki{K)p\K') J e'^^+^"^-''d^x = VoJ2^iif^)P\-f^) 

k,k' k 

we obtain Poisson brackets 

{ki{K),p\K')} = STx^GV^-Hi5^,^^. 

•^As stated above, this is sufficient for a perturbative treatment of inhomogeneities. 
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between the modes. Link integrals of the connection and fluxes then become 

= 2 ^ ~ki{K)e'''-''e''''^°/^ sin(«;j£o/2) /ki 



and 

Fjy 



r rh)/2 /.£o/2 

= 4^p"^(K)e^'^-^'e^'^^^°/'sin(Kx4/2) sin(Ki4/2)/KK«L 

where the values of indices K and L are defined such that ejxL — 1- Thus, 
{J,,„ Fj,,,} = 647r7G'V'^/ ^"'^'""'^ sin(«;,£o/2) sin(«;^£o/2) sin(fi;z,£o/2)/(«;,«;^«;i) 

= 87r7G'5/x^„(^;-^;') = 87r7G'5/V- (8) 

Using Fourier series leads to the appearance of characteristic functions Xioi^) centered at 
X of width £o. Restricted to vertices on a lattice of spacing io, this is identical to Sx,o- 
When quantized, holonomies will again become multiplication operators 

hi,v \ • • • , l^i,v, . . .) = I . . . , fii,v + 1, • • •) (9) 

and change the link labels when acting on a state. Note that we do not use a parameter 
like iiq since edges of the link size are distinguished and multiplying with holonomies for 
edges not integer multiples of lattice hnks does not preserve the lattice structure. The 
parameter //q is thus an integer, simply corresponding to the choice of representation of 
holonomies. For basic holonomies, //q = 1- 
Flux operators 

Fl,v\..., I^I,v, ■ ■ ■> = -jllllJ'lA I^I,v, ■ ■ ■) (10) 

corresponding to a lattice site orthogonal to an edge in direction I have eigenvalues pro- 
portional to yUj or, when we take a surface intersected by the vertex v, proportional to 
the average of labels of neighboring edges, 

^I,v\ ■ ■ ■,I^I,v, ■■■) = -^¥p{f^I,v+ei,, + fJ'I,v) \ ■ ■ ■,^ll,v, ■■■) ■ (11) 

In this notation, v + ej^^ denotes the vertex next to v along the edge ej^^, i-e. the endpoint 
of ej^y other than v. These labels then determine the volume eigenvalues. From 

V V 
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Figure 1: Lattice with elementary edge and surface. 



the vohime operator is defined as = Ylv 11/= 



with eigenvalues 




V 1 = 1 



(12) 



Also as in isotropic models, we can construct composite operators such as those for 
inverse triad components. The only difference is that we construct local versions of such 
operators localized at vertices, for which we use neighboring link holonomies. Eigenvalues 
of such commutators will then look similar to those in the isotropic case, except that we 
have a sum over vertex contributions where single link labels change by ±1 in the form 
Vv{fii,v+^) — Vv{fJ.i,v — ^), rather than one global contribution as in (jl)). Also, as noted before 
these functions do not contain /iq since the link size is fixed. By the same calculations used 
to derive d{p) in the isotropic case, such functions have peaks at values yU* ~ 1, or /i,, ^ j 
for an arbitrary irreducible representation, corresponding to ~ |7J^p- This is similar 
to the isotropic case, but now distinguishes different regimes according to i^p > or 
ilp < rather than inequalities for the isotropic p = Vq^^p. So also here, Vq has been 
replaced by the smaller £o- 



12 



4 Basic observations 



Having reviewed isotropic basic variables and introduced analogous ones for inhomoge- 
neous lattice states we can draw conclusions regarding their relation. Most of these effects 
are qualitatively true for any inhomogeneities, but are made explicit in lattice models 
introduced here. 

4.1 Rescaling freedom 

As noted before, coordinate background structures to define homogeneity or the choice of 
lattices in inhomogeneous states introduce parameters such as the cell volume Vq or the 
lattice spacing Iq. The parameters then occur in derived expressions, and such expressions 
change when other choices for the parameters arc made. Since the parameters are related to 
coordinates or the embedding of structures in space, it is not always guaranteed that such a 
rescaling freedom makes physical sense. Indeed, for Vq this is not the case which is one of the 
difficulties in purely homogeneous models. Ratios such as = p/p* = V^^^p/p^ which 
demarcate classical from quantum behavior, where ji^ is a characteristic scale appearing in 
the construction of operators and fi a state label, depend on Vq through p. (Factors of Vq 
do not cancel since p* is defined independently of Vq.) 

The situation is, however, different for which replaces Vq in elementary expressions 
when lattice states are used. Although Vq is still present, it always appears in combination 
with £o through N. Then, ratios as before take the form = p/p^ = (-Qp/p* which is 

£o-dependent. But unlike Vq, changing £q has physical meaning because we then change 
the scale on which we probe space by the lattice. If we choose a bigger lattice spacing, 
i.e. enlarge £o, p will have to drop to even smaller values for quantum corrections to 
become noticeable. Since both io and p are coordinate dependent, this statement about 
the relation between changing io and corresponding changes in p is invariant under rescaling 
of coordinates. This happens in such a way that the change of io implies physical properties 
that are expected. Isotropic models show similar technical features, but due to additional 
backgrounds involved they are not as physically convincing. 

These observations also indicate that the triad scale at which quantum effects become 
important in inverse powers does not only depend on a representation label chosen for the 
quantization of inverse triad operators but also on the lattice size. In fact, choosing finer 
lattices has the same effect as choosing higher representations on a lattice of unchanged 
size. The importance of such effects is determined by ratios p/p* where p is a flux value and 
is proportional to j. Technically, this ratio can be made small by choosing larger j, or by 
making the relevant fluxes take on smaller values in the same situation. In homogeneous 
settings one only has the first option, while inhomogeneous ones easily allow the second 
one by choosing smaller surfaces for integrating fiuxes. 

In fact, if the geometry is nearly isotropic, we can use the average value p — Vq^^P 
to make contact with isotorpic variables. Now, as discussed at the end of the preceding 
section, quantum corrections start to become noticeable when i^p ^ p^, or p = V^^^p = 
N'^^^IqP a; N'^^^p^ which is enhanced by a factor depending on the number of vertices. Note 
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that we did not use higher representations for holonomies in commutators, by which one 
can achieve a similar effect if j ~ iV^/^. The enhancement comes just from the fact that 
for individual links it is the local value of p rather than the total one which is relevant. 

4.2 Higher spin representations 

This observation brings us to the next point, which is the naturalness or unnaturalness, 
or even consistency, of using higher representations of holonomies to construct operators. 
As we have seen, in inhomogeneous systems one can achieve the same effect by using 
finer lattices which is certainly a legitimate way to change parameters. While this is not 
available in homogeneous models, using higher spins there may simply be seen as a way 
to mimic inhomogeneous behavior in such a setting. There is no direct relation between 
all features of homogeneous models to the full theory because of degeneracies: changing 
very different ingredients of an inhomogeneous situation can result in the same change in 
a homogeneous model. The higher spin compared to the lattice size is one such example, 
where in earlier papers only the direct relation identifying a higher representation in a 
homogeneous model with a higher representation in the full theory has been made. If this 
is the only way to relate higher spins to properties of the full theory, it certainly makes 
higher spins in homogeneous models look unnatural. The relation to the lattice spacing, 
overlooked so far, makes using higher representations much more natural in homogeneous 
models, while inhomogeneous ones can still be formulated by restricting oneself to the 
fundamental representation in composite operators. 

Along similar lines one can justify using different spins for gravitational and matter 
parts in a Hamiltonian constraint as effective means to include inhomogeneous features. 
Even if we use higher spins in homogeneous models, one may still question why one should 
use different ones in holonomies in a matter Hamiltonian compared to the gravitational 
part of the constraint. Such different spins are often helpful to bring out matter corrections 
as the dominant ones since they are most easy to derive and to work with. The relation to 
lattice spacing suggests a reason for higher spins in matter terms compared to gravitational 
ones: it simply means that for matter terms we use longer range "interactions" based on 
holonomies which extend through several vertices rather than just one basic lattice link. 
This is an option one clearly has in any lattice model, for which loop quantum gravity 
is one example. If one accepts such longer range interactions, one is led directly to the 
behavior also resulting from higher spins in a homogeneous model. 

4.3 Cosmological constant 

In addition to ratios p/p* relevant for inverse triad operators, connection components 
entering holonomies are the second variable which determines where quantum corrections 
become noticeable. In isotropic models, HqC has to be small compared to one for the classical 
constraint as a polynomial in c to be a good approximation of functions of exponentials 
exp(i/ioc/2) = exp(i//oV^^^^c/2) used for the quantization. As noted before, however, in the 
presence of a cosmological constant c can become arbitrarily large even in classical regimes 
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(while fiQ has been argued to be of order one by comparing with the lowest area eigenvalue 

m)- 

In inhomogeneous lattice models, holonomies appearing in composite operators are 
associated with links and of the form exp(iioc/2) = exp{iN~^^^c/2) where we introduced 
the isotropic c = Vf^^^c for comparison. What has to be small is thus not fiQC but N~^^^c 
which can be achieved by having a large number of vertices even if c does not become 
small. A large number of vertices is in fact what one expects if a universe grows to larger 
size because the fundamental Hamiltonian of loop quantum gravity most likely creates new 
vertices in its action |32] . The creation of new vertices is not easily incorporated in a 
regular lattice model, although it is possible, but one can view regular lattices as a coarse- 
grained, effective description of the fundamental one. For this, no explicit mechanism 
is required: for effective equations we can compute expansion coefficients locally on a 
given lattice with N vertices. In those coefficients, metric values appear depending on 
N. Parametrically, N can then be assumed to be a function of the total volume or scale 
factor whose precise behavior is to be determined by a more detailed formulation. This is 
analogous to a mean field approximation where "fast" modes are treated only effectively 
by correction terms. Here, vertices are created at each step of the action of the constraint 
operator implying Planck-size changes in the local geometry. This degree of freedom is 
thus much faster than semiclassical effects we are mainly interested in here, i.e. it happens 
on much smaller time scales (measured in terms of the global geometry). 

Then, as the total volume increases one should also adapt the lattice size and increase 
A^. This is in particular the case when one has a positive cosmological constant because 
the universe grows without bound. As the scale factor a increases and c oc ^/Xa increases, 
also A^ has to increase. For a well-behaved semiclassical state this has to happen in such 
a way that c/N^^^, the argument of holonomies, remains small which is easily achieved in 
an inhomogeneous model. The validity of a homogeneous model, on the other hand, will 
break down at a certain volume, which is one example where homogeneous approximations 
are good on small scales but not at large ones.'* 

We can draw one further conclusion related to this observation. In link holonomies 
we have the quantity N~^/^c compared to the isotropic ixqc with a constant /iQ. In inho- 
mogeneous situations, /xq does not appear at all, and there is thus no indication for the 
makeshift argument relating fiQ to the area spectrum used in the absence of a relation to 
inhomogeneity. This agrees with the fact that the full constraint operator does not contain 
the area operator. Curvature components of the constraint are rather quantized directly 
through holonomies, which is also realized here with natural loops provided by the lattice. 
Moreover, if N^^^ is not treated as a constant, we can resolve problems at large scales. This 
can be modeled in a purely isotropic situation by letting /io be scale-dependent in a form 

''On very small scales close to classical singularities, of course, inhomogeneities grow and homogeneous 
models can only give qualitative insights. Classically, one would expect better and better behavior of a 
homogeneous model as an approximation to inhomogeneous behavior the larger the universe grows. This 
is not realized in quantized models with a discrete spatial structure because more and more elementary 
blocks are needed to describe a large space. Very close to classical singularities, on the other hand, most 
indications so far show that anisotropics rather than inhomogeneities are the crucial degrees of freedom. 
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fi{p) ~ IpI^^/^ at large scales, corresponding to the number of vertices growing linearly with 
volume, oc This is exactly the behavior which has been proposed within isotropic 

models to counter large scale and other problems W7\ . While one might question a 
non-constant fiQ in isotropic models on the grounds that holonomies would have to depend 
on triad variables, the relation to the number of vertices in inhomogeneous states clearly 
provides convincing justification for it. In general, however, the behavior may differ from 
oc \pf^^ which is to be checked in detailed models implementing the subdivision. There 
are also conceptual differences between the arguments for ft{p) in isotropic models and 
those presented here: while in isotropic models the parameter arises at the operator level 
and leads to holonomies in the constraint depending on triad components, lattice models 
implement a volume dependence through states. This is closer to the full theory where it is 
the state-dependent regularization of the Hamiltonian constraint that determines how new 
edges and vertices are created. Nevertheless, it is promising and encouraging for other con- 
structions in models that independent considerations in homogeneous and inhomogeneous 
settings lead to the same qualitative conclusions. 

4.4 Relative size of quantum corrections 

Finally, we discuss the relative size of different quantum corrections arising from quantum 
Hamiltonians. The main corrections are from inverse triad operators, where the size of p is 
relevant, and holonomies where c is relevant. While p has to be sufficiently small for inverse 
triad corrections to become large, c has to be large for large holonomy corrections. For 
evaluations of effective equations [33 EHl IHH 1121 HOI HH] it is then helpful to know if there 
is any correction which dominates, or if all of them have to be taken together which would 
make the analysis more complicated. In any model, an answer to this question depends 
on which regime one is looking at, and the situation in homogeneous models is undecided. 
Irrespective of what the answer in a particular homogeneous scenario is, however, the 
situation is different in inhomogeneous models. As before, inhomogeneity implies that for 
local operators appearing in composite ones connections are integrated only over small 
links and triads only over lattice sites rather than all of space. Thus, the connection 
values relevant for holonomies and the triad values for fluxes are both significantly reduced 
compared to the values that appear in homogenous operators. As a consequence, inverse 
triad corrections will become more prominent while holonomy corrections will be less so. 
Main effects are then to be expected from inverse triad corrections of quantum geometry, 
as they have been studied several times in homogeneous models (see, e.g., jSHl IMl UHl 

EziEHiEniinn]). 

5 Further applications 

Besides these technical observations we can also draw more general conclusions. 
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5.1 Symmetry reduction: from inhomogeneity to homogeneous 
models 

Symmetric states are obtained by restricting full states to a subset of invariant connections 
relevant for a symmetric model ^U]. Such states are necessarily distributional in the less 
symmetric situation, and general operators will not leave the space of such states invariant. 
Nevertheless, one can derive all basic operators given by holonomies and fluxes of the 
symmetric model from corresponding ones in the full theory. This is sufficient to derive 
the quantum representation of models from the unique one in the full theory |H2] • Since 
many properties of loop quantizations, including those discussed here, follow already from 
basic aspects the derivation of the basic representation is a crucial step. Explicit examples 
for such constructions can be found for reducing an anisotropic model to an isotropic one 
in [HS] and for obtaining the spherically symmetric representation from the full one in ^H] ■ 
A more general discussion is given in 

We can use the constructions presented here to shed more light on the role of inhomo- 
geneities in these reductions. A general lattice state is of the form 

I,v 

with kj{x) = ioki{x) and hj^y ~ exp{ifijyki{v)/2). This corresponds to an isotropic 
connection c if ki{x) = c, or kj{x) = ioVf^^^^c = N^-^^'^c, for all / and x. The restriction 
of the state then becomes the isotropic state 

ip^{c) = exp(i/ic/2) =: (c|/i) with /x = A^"^/^ ^ fij^y . (13) 

I,v 

Following the general procedure, there is a map a taking an isotropic state of the reduced 
model to a distributional state (yu| = cr{\fi)) in the inhomogeneous setting such that 

(//| = for all I . . . . . .) (14) 

where the inner product on the right hand side is taken in the isotropic Hilbert space, 
using the restricted state | . . . ,/i7,^, . . as an isotropic state (fT^ . 

Link holonomies as multiplication operators simply reduce to multiplication operators 
on isotropic states. Fluxes for lattice sites, however, do not map isotropic states to other 
isotropic ones. This can easily be seen using {n\Fjy | . . . , /x/^t,, . . .) = It^pA* jy | • • • , fJ'i,v, ■ ■ •) 
on states 

1^,,,) := |0,. .. ,0,1,0,. ..,0) 
which have non-zero labels only on one lattice link e/^^,. We then have 

{fJ'\FLv\i'i,v) = ^l^llJ'i,v5^,,i and (/i|F7,„|?/'/,^+e^_J = 

since the flux surface and the non-trivial link do not intersect in the second case. However, 
i^\4'i,v) = {j^\4'i,v+ei) for any isotropic state (z/|. Thus, {^\Fi,v cannot be a superposition 
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of isotropic distributional states, and flux operators associated with a single link do not 
map the space of isotropic states to itself. (The above formulas show that (1| cannot be 
contained in a decomposition of {fi\Fj^y in basis states, but we can repeat the arguments 
with arbitrary values for the non-zero label in \ipi^v)-) 

Instead, one can use extended fluxes which are closer to homogeneous expressions. 
First, we extend a lattice site flux F/ „ to span through a whole plane in the lattice, leading 

This corresponds to a homogeneous flux in the box of size Vq but is still not trans- 
lationally invariant because the plane {v' : Vj = vj} is distinguished. We can make it 
homogeneous on the lattice by averaging along the direction / transversal to the plane. 
This leads to a sum over all lattice vertices in 

V 

including a factor A^~^/^ from averaging in one direction. Finally, we take the directional 

average 

I I,v 

to define the isotropic flux operator. 

Now, to find the action of this operator on a distributional state we compute 

. . . , . . .) = /ijy (z^l • • • , ^^i,v, ■■■) = 

J,v' 

where fi is defined in terms of as in (fT!?jl . This agrees with the isotropic flux operator 
defined in isotropic models, 

pailfi)) = (T{p\fi)) , (15) 

and shows in particular that p, unlike p/^^, maps an isotropic distribution to another such 
state. Thus, the isotropic representation in loop quantum cosmology follows from the 
inhomogeneous one along the lines of a symmetry reduction at the quantum level. Notice 
that this leads directly to an operator for p rather than p without explicitly introducing 
the cell size Vq. We just extend fluxes over the whole lattice in order to make them 
homogeneous, such that the combination V^^^p = p automatically arises. 

5.2 Gravitons as collective excitations 

On inhomogeneous lattice states we can look for representations of inhomogeneous excita- 
tions of geometry. Classically, such excitations are given by modes of a metric perturbation 
lab = Qab + hab ou a backgrouud (jab- In a perturbative quantization on that background, 
hab would thus be the basic field to be quantized, leading to the expectation of one-particle 
states called gravitons corresponding to transverse-traceless modes of hab- 
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In loop quantum gravity, the situation is different because its basic formulation is 
background independent. Having formulated here a lattice description as a sector which 
re-introduces a metric background, we can see how such metric perturbations are realized. 
Simplifications due to Abelianization used here also occur in constructions of explicit states 
in linearized gravity [HZj. As in the full theory, basic excitations of the loop quantization 
are given by dynamical moves which change the labels /i/^^ locally at a vertex to fij^^ ± 1. 
This changes the local labels, corresponding to inhomogeneous modes, but also the total 
volume corresponding to the background volume of Qab- A quantum analog of the mode 
can be defined as the difference 

Spi,v = Pl,v -P= \l^l^I,v - (K({/i/,4)/Ar)2/3 

between local labels and the average volume which captures properties of the classical 
modes even at the level of effective dynamics pT] . 

The quantum picture of these excitation is then very different from the classical one: 
6pi^v arises non-locally since we have to sum over the whole lattice to subtract off the 
background. Basic dynamical moves of the quantum theory change the local quantum ex- 
citation as well as the background, leading to a mode change also in a non-local way. In this 
picture, metric modes such as gravitons arise as collective excitations out of the underlying 
lattice formulation: many basic quantum modes combine to give classical excitations. 

5.3 Quantum gravity corrections for large universes 

Quantum expressions differing from their classical analogs give rise to corrections to clas- 
sical equations. This usually occurs on small length or large curvature scales as they are 
realized in the very early universe. But subdivision in a lattice also makes relevant lo- 
cal length scales smaller, and indeed we already noticed that inverse power corrections as 
functions of p\v) become dominant over holonomy corrections because the flux values are 
reduced on lattice links if the lattice becomes finer. This opens the possibility for small 
corrections on single lattice sites to add up to sizeable corrections on the whole lattice, 
which can influence the evolution of modes [21] or even of the whole universe. 

6 Conclusions 

Models of loop quantum gravity are being investigated actively regarding their phenomeno- 
logical properties, and perturbative inhomogeneities are currently being included. This 
brings us closer to reliable computations of potentially observable properties such as those 
of structure formation. It is then important to check all intrinsic details of such models 
and see how faithfully they incorporate features of the full theory. As we discussed, quali- 
tative effects are realized in homogeneous as well as inhomogeneous lattice models in the 
same way. We introduced inhomogeneous lattice constructions in a way which allows for 
a relation to a homogeneous background. The relation to isotropic models is then clear. 
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which provides a new step toward relating isotropic models to the full theory. Although 
the lattice construction is analogous to that of isotropic models, quantitative aspects can 
change which has a bearing on which ranges of parameters one considers as natural or 
unnatural. It also plays a role for which correction terms will be dominant in different 
regimes which is the most important aspect for phenomenological investigations. 

The lattices used here are not intended as fundamental descriptions since they are not 
background independent. They are rather to be thought of as effective lattice models 
which result if some degrees of freedom of the full theory are interpreted as providing a 
background for other degrees of freedom. This simplifies the more irregular fundamental 
behavior. It is this feature which helps to illustrate and clarify some puzzling properties 
realized in isotropic models. There are also further applications such as explicit construc- 
tions of quantum field theories on a quantum spacetime as suggested in jHHl 120], more 
general models for field propagation with quantum corrections based on 1^3] , or cosmolog- 
ical perturbations j2I]. Such lattice models are more accessible for a relation to the full 
theory and can thus be used to suggest relevant properties which one wants to ensure in 
future constructions and specifications of the full theory at a basic level. 

For now, the construction already presents several applications. They all rely on the 
re-introduction of a background which is necessary to interpret computational results. We 
focused here on kinematical aspects, but an inclusion of characteristic dynamical ones is in 
progress. For a detailed description of all effects expected in fully inhomogeneous settings, 
remaining difficulties are an explicit treatment of subdivision and the inclusion of non- 
Abelian properties. But this is not required explicitly for effective perturbations where the 
lattice description discussed here is sufficient. It allows one to derive perturbation equations 
for metric modes around the background, to be used for instance for cosmological structure 
formation, or to understand the emergence of graviton states conceptually. 

Note added: After this paper had been submitted to the journal General Relativity 
and Gravitation, the preprint |gr-qc/0607100 by Kristina Giesel and Thomas Thiemann 
was posted. Although the conceptual setup is quite different from the approach followed 
here, it is clear that there is a convergence of different ideas pursued recently in loop 
quantum gravity. Especially in combination with the results of Sec. 15.11 there is now close 
contact between the full theory (understood as a general framework to set up background 
independent quantum kinematics and dynamics without symmetry assumptions) and loop 
quantum cosmology. 
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